Temperature- and magnetic-field-dependent thermal conductivity of 
pure and Zn-doped Bi2Sr2CaCu20s+<5 single crystals 
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The thermal conductivity k of Bi2Sr2CaCu2 08+* is measured in pure and Zn-doped crystals as a 
function of temperature and magnetic field. The in-plane resistivity is also measured on the identical 
samples. Using these data, we make a crude estimate of the impurity-scattering rate F of the pure 
and the Zn-doped crystals. Our measurement show that the "plateau" in the k(H) profile is not 
observed in the majority of our Bi-2212 crystals, including one of the cleanest crystals available to 
date. The estimated values of V for the pure and Zn-doped samples allow us to compare the k(H) 
data with the existing theories of the quasiparticle heat transport in d-wave superconductors under 
magnetic field. Our analysis indicates that a proper inclusion of the quasiparticle- vortex scattering, 
which is expected to play the key role in the peculiar behavior of the k(H), is important for a 
quantitative understanding of the QP heat transport in the presence of the vortices. 

PACS numbers: 74.25.Fy, 74.62.Dh, 74.72.Hs 



Thermal conductivity k of a superconductor is one of 
the few probes which allow us to investigate the quasi- 
particle (QP) density and its scattering rate in the su- 
perconducting (SC) state. It is now believed that the 
SC state of the high-T c cuprates is primarily d x 2_ y 2 [Q, 
where it has been found that the magnetic field induces 
extended QPs whose population increases as \/H~ 
Also, the QP scattering rate in the cuprates has been 
found to drop very rapidly below T c [Q, which causes 
a pronounced peak in the temperature dependence of k 
In 1997, Krishana et al. reported an intriguing re- 
sult from their measurement on the magnetic-field de- 
pendence of k in Bi2Sr2CaCri208+5 (Bi-2212) at tem- 
peratures below 20 K Q. They observed, in the profile 
of k(H), a sharp break at a "transition field" Hk and a 
subsequent plateau region where n does not change with 
magnetic field. Krishana et al. proposed an interpreta- 
tion that Hk marks a phase transition from the d x 2_ y 2 
state to a fully-gapped d x 2_ y 2 +id xy (or d x 2_ y 2 -Ms) state 
and thus in the plateau region there are little thermally- 
excited QPs which contribute to the heat transport. This 
interpretation appears to be fundamentally related to the 
high-T c mechanism and therefore attracted much atten- 
tion both from theorists [f7j-[l4| and from experimentalists 

mi- 

An independent test of this unusual behavior of k has 
been reported by Aubin et al. Q ; although the plateau- 
like feature was essentially reproduced in their measure- 
ment, when the data were taken with field sweep up and 
down, Aubin et al. observed a rather large jump in k 
upon field reversal and consequently the k(H ) profile had 
a pronounced hysteresis. The fact that the k value in the 
"plateau" depends on the history of the applied magnetic 
field casts a serious doubt on the Krishana et aVs inter- 



pretation. Moreover, Aubin et al. reported that an in- 
crease in k with magnetic field was observed at subkelvin 
temperatures, which strongly suggests the presence of a 
finite density of QPs at low temperatures and thus is in- 
compatible with the fully-gapped state p6| . Although 
these newer results suggest that a novel phase transition 
in the gap symmetry is not likely to be taking place, the 
plateau in the n(H) profile and the sensitiveness of k on 
the magnetic-field history are still to be understood. One 
interesting information one can draw from these exper- 
iments is that phonons are not scattered by vortices in 
cuprate superconductors ||[l5 19 1. 

Motivated by these experiments, there appeared sev- 
eral theories that try to capture the essential physics of 
the QP heat transport in the d-wave superconducting 
state. It has become rather clear that a proper inclu- 
sion of the QP-vortex scattering rate is necessary for ex- 
plaining the observed magnetic-field dependence; how- 
ever, there has been no consensus yet as to how the QP- 
vortex scattering is taken into account. To improve our 
understanding of the QP heat transport under magnetic 
field, quantitative examinations of various theories in the 
light of actual data are indispensable. For Bi-2212, how- 
ever, previously published data from various groups do 
not provide enough information; for example, the impu- 
rity scattering rate, which is the most important param- 
eter that controls the thermal conductivity behavior, are 
not known for clean Bi-2212 crystals. (Rather surpris- 
ingly, no electrical resistivity data have been supplied for 
samples which were used in recent studies of n(H ) profile 
of Bi-2212 §l|g@.) 

In this paper, we present results of our measurements 
of k(T) and k(H ) of well-characterized Bi-2212 crystals, 
together with their in- plane resistivity (p a b) data. To 
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look for the effect of changing impurity-scattering rate, 
we measured both pure and 0.6%-Zn-doped crystals. The 
crystals used here are single domained (without any mo- 
saic structure nor grain boundaries) and have very good 
morphology, which were confirmed by x-ray analyses and 
the polarized-light optical microscope analysis. In the 
data presented here, neither the pure sample nor the Zn- 
doped sample show any plateaus in the k(H) profile in 
the temperature and the magnetic-field regime where the 
plateaus have been reported. In fact, we have found that 
the plateau in the k(H) profile is not a very reproducible 
feature (we have thus far observed the plateau-like fea- 
ture in only 2 samples out of more than 30 samples mea- 
sured) and we have not yet conclusively sorted out what 
determines the occurrence of the plateau |2Cf] ; we there- 
fore decided to show only the data that are representa- 
tive of the majority of the samples. Using all the data of 
Pab{T), k(T), and n{H) 1 we try to estimate the electronic 
thermal conductivity K e (T) and make a rough estimate of 
the impurity scattering rate T for the pure and Zn-doped 
samples. Our data offer a starting point for the quan- 
titative understanding of the QP heat transport in the 
e?-wave superconductors in the magnetic fields, where an 
interplay between the QP-vortex scattering and the QP- 
impurity scattering apparently plays an important role. 

The single crystals of Bi-2212 are grown with 
a floating-zone method and are carefully annealed 
and quenched to obtain a uniform oxygen content 
pjj . Both the pure and the Zn-doped crystals 
[Bi2Sr2Ca(Cui_ z Znz)208+,5] are tuned to the optimum 
doping by annealing at 750 °C for 48 hours, after 
which the transition width of about 1.5 K (measured 
by the dc magnetic susceptibility measurements) is 
achieved; this indicates a high homogeneity of these 
crystals. The Zn-doped crystal contains (0.6±0.1)% of 
Zn (namely z=0.006±0.001), which was determined by 
the inductively-coupled plasma spectrometry. The zero- 
resistance T c 's are 92.4 K for the pure crystals and 84.5 K 
for the Zn-doped crystal. We first measure the tempera- 
ture dependence of p a b using a standard four-probe tech- 
nique, and then measure the thermal conductivity n of 
the same sample. The temperature dependence of k from 
1.6 to 160 K are measured in zero field using calibrated 
AuFe-Chromel thermocouples. The precise magnetic- 
field dependence of k is measured with a steady-state 
technique using a small home-made thin-film heater and 
microchip Cernox thermometers. The bottom end of the 
sample is anchored to a copper block whose temperature 
is carefully controlled within 0.01% stability and accu- 
racy. Since we need to measure the change in n with a 
very high accuracy, the small magnetic-field dependence 
of the thermometers were carefully calibrated beforehand 
using a SrTiC>3 capacitance sensor and a high-resolution 
capacitance bridge. As a check for the measurement sys- 
tem and the calibrations, we first measured the thermal 
conductivity of nylon and confirmed that the reading is 



indeed magnetic-field independent within 0.1% accuracy. 

Figure 1(a) shows the p a b(T) data of the pure and Zn- 
doped samples. If we define the residual resistivity po by 
the extrapolation of the T-linear resistivity to K, we 
obtain slightly negative po for the pure sample; this is 
always the case for the cleanest Bi-2212 crystals grown 
in our laboratory. As is expected, the Zn-doped sample 
gives larger po, which is about 10 pttcm. We note that 
the uncertainty in the absolute magnitude of p a b and k in 
our measurements are less than ±5% in this work. This 
is achieved by determining the sample thickness, which is 
usually the main source of the uncertainty, by measuring 
the weight of the sample with 0.1 pg resolution. We used 
relatively long samples here (the length of the pure and 
the Zn-doped samples were 4.5 mm and 6.5 mm, respec- 
tively), with which the errors in estimating the voltage 
contact separation and the thermocouples separation are 
less than ±5%. It should also be noted that a sophisti- 
cated technique to make the current contacts uniformly 
on the side faces of the crystals are crucial for reliably 
measuring p a b of Bi-based cuprates. 

Figure 1(b) shows the temperature dependence of k of 
the two samples in zero field. To our knowledge, the size 
of the peak in k(T) of this pure sample is the largest ever 
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FIG. 1. T dependence of (a) p a b and (b) n of pure and 
0.6%-Zn-doped Bi-2212 crystals in zero field. The thin solid 
lines are T-linear fits to the p a b(T) data. 
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reported for Bi-2212 (the enhancement from the mini- 
mum near T c to the peak is 16%). This implies that 
the pure crystal reported here is among the cleanest Bi- 
2212 crystals available to date. The Zn-doped sample 
shows not only a smaller magnitude of k but also a signif- 
icantly suppressed peak in k(T) below T c ; this is caused 
by a larger impurity scattering rate (which limits the en- 
hancement of the QP mean free length below T c ) and 
is consistent with the Zn-doping effect in YBa2Cu307_^ 
(YBCO) (2§§|. 

Figure 2 shows k(H) profiles of the two samples at 7.5 
and 12.5 K. The data are normalized by the zero-field 
value of k. Neither of these samples show any plateau- 
like feature below 6 T. (The k(H ) data of the Zn-doped 
sample is almost flat within our sensitivity, but we do 
not call it a plateau.) Note that Ref. || reported the 
plateaus to be observed above 0.9 T at 7.5 K and above 
2.6 T at 12.5 K, which are well within the range of our 
experiment. The data in Fig. 2 were taken in the field- 
cooled (FC) procedure, in which the sample is cooled in 
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FIG. 2. Normalized k(H) profiles of the pure and Zn-doped 
samples at (a) 7.5 K and (b) 12.5 K. The solid lines are the 
fits to Eq. (1) with: (a) K ph =2.60 W/Km, K e =0.10 W/Km, 
and p=0.88 for the 7.5 K data, and (b) n ph =2.91 W/Km, 
K e =0.20 W/Km, and p=0.53 for the 12.5 K data. 



the presence of a magnetic field, and therefore the mag- 
netic induction in the sample remains homogeneous. The 
FC data are expected to give information which is free 
from complications due to vortex pinning, while the zero- 
field-cooled (ZFC) data (for which the sample is cooled 
in zero field and the magnetic field is swept at a constant 
temperature) are subject to such complications. We em- 
phasize that the FC data should be examined to look for 
the intrinsic effect of vortices on the QP transport. All 
the previously published data on the plateau P, p^| , ^6| , ^8[ 
are, however, ZFC data pi| . 

It has been proposed that the magnetic-field depen- 
dence of k of the high-T c cuprates are described by 



k{H,T) = 



l+p(T)\H\ 



(1) 



where K e is the electronic part of k in zero field and K p h is 
the phonon part [p6| . Equation (1) was proposed first by 
Yu et al. [^7| and later by Ong and co-workers J5|, ^9| , p5]| . 
This expression utilizes the finding that the phonon ther- 
mal conductivity of the cuprates is independent of the 
magnetic field. The parameter p(T) is proportional to 
the zero-field value of the QP mean free path. We found 
that the k(H) data of our clean samples are reasonably 
well described by Eq. (1). The solid lines in Fig. 2 (a) 
and (b) are fits of the data to Eq. (1). The fitting pa- 
rameters suggest that the phononic contribution to n is 
as large as 96% and 93% at 7.5 and 12.5 K, respectively; 
this is essentially a reflection of the fact that the changes 
in k with the magnetic field are very small at these tem- 
peratures. Ong et al. reported |l9| that p{T) is about 
2.3 T- 1 for underdoped YBCO at 7.5 K, while it is 0.88 
T _1 for our pure Bi-2212 at 7.5 K. This is an indication 
that the QP transport is dirtier in Bi-2212 compared to 
YBCO; we will elaborate on these p(T) values later. The 
magnetic field dependence of n of the Zn-doped sample 
is too small to make a reliable fit to Eq. (1), but one can 
infer that the ratio of the phononic contribution in the 
Zn-doped sample is even larger than in the pure sample. 

In the normal state, the electronic thermal conduc- 
tivity n e and the electrical conductivity a are related by 
K e /T = La, where L is called the Lorenz number. In sim- 
ple metals, L is usually constant at high temperatures 
(Wiedemann-Franz law) and the free-electron model 
gives L = 2.44xl0~ 8 WQ/K 2 . When the electron- 
electron correlation becomes strong, L becomes smaller 
than the free electron value Ln; for YBCO, L has been es- 
timated to be around 1.0xl0~ 8 WQ/K 2 near T c 
Using this value of L and the p a b data of our crystals, 
we can roughly estimate n e and k p h above T c . Such es- 
timation for the pure sample gives K e and n p h values of 
about 1.6 and 3.6 W/Km, respectively, at 120 K. (120K 
is the lower bound for the temperature range where the 
effect of the superconducting fluctuations on p a b is negli- 
gible.) It has already been established that K p h does not 
change much just below T c and that it is mostly K e that 
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causes the peak J5|,|29| . We can therefore infer (from the 
K e value at 120 K and the total k at the peak) that n e 
at the peak is about 1.7 times larger than in the normal 
state. In the same manner, we can estimate for the Zn- 
doped sample K e and k p h to be about 1.0 and 2.9 W/Km 
at 120 K, respectively, and the enhancement of K e at the 
peak is inferred to be a factor of 1.2. This analysis indi- 
cates that the effect of 0.6%-Zn doping is weaker for K p h 
(which is decreased by 19% at 120 K upon the 0.6%-Zn 
doping) than for K e (which is decreased by 38% at 120 
K). 

Comparison of the inferred behaviors of n e below T c 
with the theoretical calculations of Ref. 0] gives us the 
idea of the magnitude of the impurity-scattering rate T in 
our samples. As is discussed above, n e can be inferred to 
be about 1.7 times enhanced at the peak; comparison of 
this enhancement factor with the calculations for various 
r/T c values suggests that T/T c of our pure sample is 
about 0.05. The position of the peak (at T/T c ~ 0.7) is 
also consistent with the theoretical calculation for T/T c ~ 
0.05. Similarly, the inferred enhancement factor of K e {T) 
of the Zn-doped sample suggests T/T c to be about 0.2. 
Although these estimates are very crude, the estimated 
values of T/T c of our Bi-2212 samples imply that the 
scattering rate increases by T /T c ~ 0.25 per 1% of Zn, 
which is of the same order of magnitude as that for YBCO 
HQ. The estimated T/T c of our pure Bi-2212 is still 
notably larger than that for pure YBCO, for which T/T c 
has been estimated to be ^0.01 [^2],[50) . This is essentially 
a reflection of the fact that the peak in k(T) of Bi-2212 is 
much smaller compared to that of YBCO, and is probably 
caused by the intrinsic disorder of the crystalline lattice 
of Bi-2212 (the modulation structure along the b axis). 

The above mentioned difference in T/T c between our 
pure Bi-2212 (T/T c ~ 0.05) and pure YBCO (T/T c ~ 
0.01) indicates that the QP mean free path in zero field, 
Iq, is roughly 5 times longer in pure YBCO compared to 
that in pure Bi-2212. This observation yields a useful 
information on the QP scattering cross section of the 
vortices, <Jt r , in Bi-2212. It was discussed in Refs. Jl]]] 
and [^BTl that parameter p(T) appearing in Eq. (1) can 
be expressed as p(T) = loatr/^o- (Note that a t r is the 
scattering cross section in two dimensions.) As is already 
discussed, the p(T) value at 7.5 K is 0.88 T _1 for our 
pure Bi-2212, while it is 2.3 T _1 for YBCO; given the 
indication that lo is roughly 5 times longer in YBCO, 
the ratio of the p(T) values of the two systems suggests 
that atr should be approximately a factor of 2 larger in 
Bi-2212. Since atr has been estimated to be 9 nm for 
YBCO [0, we obtain a tr ~ 18 nm for Bi-2212. The 
origin of this difference in a tr between the two systems 
might be related to the difference in the structure of the 
vortex lines and should be a subject of future studies. In 
any case, the estimate of atr ~ 18 nm gives Zo ^ 0.1 /im 
for pure Bi-2212 at 7.5 K. 

Now let us discuss the observed magnetic-field depen- 



dence of k in conjunction with the estimated values of 
r/T c for each sample. Kiibert and Hirschfeld (KH) cal- 
culated the magnetic-field dependence of K e that comes 
from the QP's Doppler shift around the vortex 
Good agreements between the KH theory and experi- 
ments have been reported for very low temperatures, 
where an increase of k with magnetic field has been 
observed |l^,^,^]. The numerical calculation by KH 
show that k(H) is already an increasing function of H 
at T=0.2T C for a dirty case, r/T c =0.1; this is clearly in 
disagreement with our data and indicates the necessity 
of an inclusion of the QP-vortex scattering into the cal- 
culation. 



The theory proposed by Franz 13 tries to incorporate 
the effect of QP-vortex scattering; heuristically, Franz 
supposed that the QP-impurity scattering and the QP- 
vortex scattering are separable and additive. In this the- 
ory, the zero-field scattering rate ao (which is a sum 
of the impurity scattering rate and the inelastic scat- 
tering rate) is directly related to the total change in 
K with the magnetic field, Ak, which is expressed as 
Ak = (2.58T/a Q - l)n o for a < T, and as Ak = 
(2.15T/ao) 2 KQQ for ao > T. In these expressions, koq 
is the universal thermal conductivity |50 33 1. Assuming 
that the inelastic scattering is negligible (ao w T) at 7.5 
K, we can estimate Ak for the pure and Zn-doped sam- 
ples using these equations. The results are Ak sw 0.4 
W/Km for the pure sample and Ak 0.1 W/Km for 
the 0.6%-Zn-doped sample. (In the calculation, we used 
Kqo/T = 0.015 W/K 2 m which is reported for Bi-2212 
p2[.) When we turn to our data, the values of Ak ac- 
tually measured in our experiments are much smaller; at 
7.5 K, Ak up to 6 T is 0.086 and 0.008 W/Km for the 
pure and the Zn-doped samples, respectively. Clearly, the 
theory overestimates Ak and the overestimation is par- 
ticularly large for the Zn-doped sample. This comparison 
implies that the rather simple treatment of the separa- 
ble QP-impurity scattering and the QP-vortex scattering 
is not accurate enough, particularly when the impurity 
scattering rate becomes comparable to the vortex scat- 
tering rate. 

Vekhter and Houghton (VH) have recently pro- 
posed a theory which explicitly considers the inter- 
play between the vortex-lattice scattering and the im- 
purity scattering. Unfortunately, numerical calculations 
for the magnetic- field dependence of k at intermediate 
temperatures are available only for a very clean case 
(r/T c =0.006) [O], and we cannot make a direct com- 
parison to our data. Qualitatively, however, the theory 
predicts that a steep drop in k(H) at low fields is ex- 
pected only when T is larger than T. The estimated T/T c 
values of our samples suggest that the condition T > T 
is satisfied in the pure sample, but is not satisfied in the 
0.6%-Zn-doped sample in the temperature range of Fig. 
2. 

In summary, the temperature and magnetic-field de- 
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pendences of the thermal conductivity k are measured 
in pure and Zn-doped Bi-2212 crystals that are well- 
characterized. The p a b{T) data taken on the identi- 
cal samples are used for extracting the electronic ther- 
mal conductivity n e above T c . The temperature and 
magnetic-field dependences of K clearly reflect the dif- 
ference in the impurity-scattering rate T in the crystals. 
It is found that in the majority of our Bi-2212 crystals 
(including one of the cleanest crystals available to date) 
the "plateau" in the k(H ) profile is not observed and the 
k(H) profile around 10 K are reasonably well described 
by Eq. (1), which was originally proposed for YBCO. 
We estimate T of our crystals by comparing the behav- 
iors of K e below T c with the calculation by Hirschfeld and 
Putikka [^2| . The estimated values of T for the pure and 
Zn-doped samples allow us to compare the k(H) data 
with the existing theories of the QP heat transport in 
d-wave superconductors under magnetic field. 

We thank H. Aubin, K. Behnia, M. Franz, T. 
Hanaguri, A. Maeda, Y. Matsuda, and N. P. Ong for 
fruitful discussions. 
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